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Abstract. 

Positive and negative pulsar breaking indices suggest that some fraction of the pulsar spindown torque undergoes 
a cyclic evolution. The observed strong correlation of 'anomalous' breaking indices with pulsar age implies that the 
characteristic periodicity timescale is in the range 100 to 10,000 years depending on the fraction of the spindown 
torque that undergoes cyclic evolution, 1 to 100% respectively. We argue that the longest variability timescale is 
consistent with a neutron star magnetic cycle similar to the solar cycle. 
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1. Observational evidence 

A pulsar spins down due to the torque on the neutron 
star crust generated by the electric current flowing in 
its magnetosphere. In the simplified picture of a steady- 
state axisymmetric force-free ideal MHD magnetosphere, 
Contopoulos, Kazanas & Fendt 1999 (hereafter CKF) first 
showed that the distribution of the magnetospheric elec- 
tric current I can be determined as an eigenvalue of the 
problem, if one makes the natural assumption that the 
magnetosphere is smooth and continuous on the light 
cylinder (defined as the distance rj, = c/fi from the ro- 
tation axis, where f2 is the pulsar angular velocity). The 
unique electric current distribution thus obtained yields a 
unique pulsar spindown torque, and thus a unique pulsar 
spindown rate f2. This is of the same order as the value 
obtained for simple electromagnetic vacuum dipole radia- 
tion, namely 
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Here, _B„, r*, and M* are the neutron star polar mag- 
netic field, radius and mass respectively; and / is a nu- 
merical factor of order unity. It turns out that this re- 
sult is also valid in the general non-axisymmetric case 
(Bogovalov 1998; Spitkovsky 2006; Contopoulos 2007a), 
and therefore, one may use Eq. Q] to obtain an estimate 
of the polar value of the neutron star magnetic field from 
measurements of f2 and Q. 

Strictly speaking, however, the results of CKF and sub- 
sequent related work are only valid in steady-state. Thus, 
as the neutron star spins down and the light cylinder 
moves to larger and larger distances, one needs to take 
into account the evolution of the pulsar magnetosphere. 



The first thing one may assume is that the magnetosphere 
evolves through a sequence of steady-state equilibria of the 
CKF type, i.e. that it manages to readjust itself so that 
at all times, the region of closed lines extends all the way 
to the light cylinder, and the last open magnetic field line 
extends to infinite distances without reconnecting accross 
the equator. In addition, one may assume that does 
not evolve with pulsar age. Unfortunately, the situation is 
more complicated than that, since Eq. [1] yields a braking 
index value 



= 3 



(2) 



and most known measurements of n differ from that value. 
In fact, there exist today ~ 400 pulsars in the ATNF cat- 
alogue (Manchester et al. 2005) with measured values of 
n in the range from -10 6 to +10 6 . Although all but six of 
these values are characterized in the literature as 'anoma- 
lous', one thing is certain: Eq.[2]cannot be right in general. 

In Figure 1 we plot ±log|n (with ± according to 
whether n > or n < respectively) as a function of 
the characteristic spindown time r = — 0/(2f2) in years. 
One may argue (Alice Harding, personal communication) 
that in young pulsars (r < 10 5 ), braking index measure- 
ments may be 'corrupted' by neutron star glitches. On the 
other hand, in older pulsars (t > 10 5 ) where glitches are 
not as important, one finds the correlation 



± log |n| ~ ±(logr — 3) 



(3) 



On the other hand, in a P — P diagram (P = 27r/0 is the 
pulsar period) containing all the cases with observed val- 
ues of n, there is no obvious correlation between P and P 
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(Figure 2). In other words, the pulsars with measured val- 
ues of n have been taken randomly from the total number 
of known pulsars without any obvious observational bias, 
and therefore, they represent pulsars at various stages of 
their evolution. Beskin, Biryukov & Karpov (2006) argued 
that the ± symmetry of Figure 1 implies that some frac- 
tion of the pulsar angular velocity undergoes a cyclic evo- 
lution. We suggest that it is more constructive to consider 
a cyclic component in the evolution of the pulsar spindown 
torque (which is, after all, the source of the angular veloc- 
ity evolution). Interestingly enough, as we will see next, 
the data also yields some information on the periodicity 
timescale. 

There are several ways to reconcile eqs. [1] and [31 (A) 
Assume that the approximation of a sequence of CKF-type 
steady-state magnetospheric equilibria holds, and that the 
neutron star magnetic field undergoes a cyclic evolution; 
(B) Relax the assumptions of the CKF analysis and as- 
sume a variable magnetospheric structure that would yield 
a cyclic evolution of the factor /; (C) Relax the assump- 
tion of constant neutron star moment of inertia. Only (A) 
and (B) refer to the spindown torque itself. In any case, 
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with f a , B* characteristic values of the spindown param- 
eter / and the polar magnetic field £?* respectively, and 
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characterizing the cyclic variation of either / or B,. Here, 
Tcyde is the characteristic period of the cyclic spindown 
evolution in years; a is the fraction of the spindown torque 
that varies periodically (0 < a < 100%); and is a ran- 
dom initial (at pulsar birth) phase angle. Note that, in 
order for the star to continuously spin down we must have 
F > at all times. Therefore, 



n = 3 + ~ 3 + 2ira 
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For old pulsars with T>T cyc i e , Eq. [6] can equivalently be 
written as 

± log ~ ±(logr - logTcyde + 1 + log a) . (7) 

Comparing eqs. [3] and one obtains the following approx- 
imate relation between the characteristic period and the 
fraction of the pulsar spindown which varies in a cyclic 
way, 



logr, 



cycle 



4 + log a 



(8) 



We plot in Figure 3 what Eq.[S]yields for the ~ 400 pulsars 
with measured values of n, assuming a = 100%. Note that 
the fit is independent of f a and -B* - As Beskin, Biryukov 
& Karpov (2006) suggested, the minimum characteristic 
period cannot be smaller than the pulsar observation pe- 
riod of 40 years, and this together with Eq.[5]yields a range 
100 < T cyc i e < 10, 000 for 1% < a < 100% respectively. 



2. Cyclic magnetospheric evolution 

Several physical models that address the issue of cyclic 
variation of the pulsar spindown have been proposed in 
the literature, ranging from neutron star interior "wob- 
bling" on a timescale of a few years (e.g. Kundt 1988), 
to magnetospheric variability (e.g. Contopoulos 2005). In 
the present work, we would like to focus on our simplest 
(one-parameter) fit of the anomalous braking index data, 
namely the one with r cyc / e ~ 10, 000 and a « 100%. 

F(t) becoming zero periodically is not compatible with 
a cyclic evolution of the neutron star moment of iner- 
tia (case C above). On the other hand such a scenario 
is compatible with a cyclic evolution of the neutron star 
magnetic field similar to the solar cycle (case A above). 
Interestingly enough, the ten thousand year timescale 
that we obtain is comparable to the neutron star cooling 
timescale (e.g. Blandford, Applegate & Hernquist 1983). 
It is conceivable that some sort of dynamo mechanism in 
the neutron star interior, may support a cyclic evolution 
with 



5* = 5*o sin(27r 
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5,000 years 



(9) 



(Eqs. |3J [5]). Note that this scenario does not require 
magnetic field decay (at least over timescales shorter 
than about 10 11 years), in agreement with the analy- 
sis of the P — P diagram presented in Contopoulos & 
Spitkovsky (2006). 

We also tried to seek variants of the CKF solution (case 
B above) that would yield values of / very different from 
unity. In fact, what we need is a physical mechanism that 
will periodically turn off the neutron star magnetospheric 
spindown. In a series of papers (Contopoulos 2007b, c), 
we relaxed the assumption of ideal MHD in the equato- 
rial region of the pulsar magnetosphere beyond the light 
cylinder. This is the region where the magnetospheric re- 
turn current flows, and several authors before us sug- 
gested that this may be the region of electromagnetic 
energy dissipation that would result in particle accel- 
eration (e.g. Coroniti 1990; Michel 1994; Lyubarsky & 
Kirk 2001; Kirk & Skjaeraasen 2003; Romanova, Chulsky 
& Lovelace 2005). As we argued in Contopoulos 2007c, 
one cannot study equatorial reconnection without taking 
into account the global topology of the poloidal magnetic 
field. The details of equatorial reconnection remain (yet) 
unknown. However, it is easy to realize that, when equa- 
torial reconnection is present, magnetic field lines that 
cross the light cylinder and would have extended to infin- 
ity in CKF, now continuously reconnect across the equa- 
tor. As a result, the equatorial condition for the magnetic 
flux function ty(r;z) (defined as the magnetic flux cross- 
ing a circle of cylindrical radius r at height z around the 
axis of rotation) differs from that in CKF. In particular, 
^(r > tl;z = 0) is not constant but decreases with dis- 
tance. We assume for simplicity that 



*(r > r L ; z = 0) = *(r = r L ; z = 0)(r/r L )~ 



(10) 
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where, e is a parameter that characterizes the effect of dis- 
sipation (e = corresponds to the ideal MHD case studied 
in CKF, whereas e = 1 corresponds to a magnetosphere 
with maximum equatorial dissipation) . Equation [TU] is a 
new (to our knowledge) equatorial boundary condition be- 
yond the light cylinder, and one may thus implement the 
same procedure as described in CKF to solve the pul- 
sar equation (Scharlemann & Wagoner 1973), and thus 
obtain the magnetosphcric structure and electric current 
distribution /('J) for various values of < e < 1, as 
seen in Figures 4-6. For each such electric current distri- 
bution the total electromagnetic spindown torque acting 
on the neutron-star crust is proportional to the integral 
/ 7(*)d* (e.g. Michel 1991). Note that when e = 0.4 
(Fig. 5) there is no equatorial return current sheet (the 
return current is distributed along the magnetic field lines 
that cross the light cylinder), whereas when e = 1 (Fig. 6) 
/ /(5')d^' ~ 0, i.e. the total neutron-star spindown torque 
is close to zero. In Figure 7 we plot the value of the spin- 
down torque parameter / as a function of our dissipa- 
tion parameter e. One sees that, as we introduce more 
and more dissipation in the equatorial region, the mag- 
netosphere evolves to a configuration with less and less 
electromagnetic torque acting on the central neutron star. 
Obviously, a cyclic evolution of the physical mechanism 
that allows or inhibits equatorial dissipation in the pul- 
sar magnetosphere (e.g. variability in the supply of charge 
carriers from the neutron star surface that may be due 
to a periodic stellar wind) would yield a cyclic evolution 
of the magnetosphcric torque. Unfortunately, magneto- 
spheric solutions for 0.4 < e < 1 (/ < 0.5, or equiva- 
lently F < 80%) contain regions with I(^f) < 0, where 
electromagnetic energy is flowing from the magnetosphere 
onto the star, and therefore, such solutions are probably 
unphysical. 

We conclude that a cyclic component in the evolution 
of the magnetosphcric spindown torque may account for 
the measured large positive and negative anomalous brak- 
ing index values. If we are willing to consider a 100% cyclic 
evolution, this can only be due to a neutron star magnetic 
cycle similar to the solar cycle. In that case, the evolution 
timescale would be on the order of 10,000 years. 
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Fig. 1. Braking index as a function of characteristic spin- 
down time. We plot here ± log \n\ (± according to whether 
n > or n < 0) vs. logr, where r = — £1/2(1 in years. 
Note that \n\ > 1 everywhere. At large r > 10 5 , the dia- 
gram may be fit by the simple linear relation =b log |n | ~ 
±(logr-3). 
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Fig. 3. Fit of the distribution shown in Figure 1 assum- 
ing a 100% cyclic evolution of the pulsar spindown torque 
(Eqs.H&El). T cycie ~ 10,000 (in years). The fit is accept- 
able, even for young pulsars (t < 10 5 ) where some of the 
dispersion in the measurements of n is due to neutron star 
glitches. 
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Fig. 2. P—P diagram for the ~ 400 pulsars with measured 
braking index values shown in Figure 1. Here, P, P are 
the pulsar period and period derivative respectively. The 
distribution is that of a standard sparse P — P diagram 
without any obvious observational bias. 
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Fig. 4. On the left, we plot the magnetic field structure 
in the case of no magnetospheric reconnection (CKF). 
Distances are normalized to the light cylinder distance 
r^. On the right, we plot the corresponding electric cur- 
rent distribution / = along the field lines that cross 
the light cylinder. The magnetic flux is normalized to the 
canonical value ^ = irB*r 3 jr^. The electric current is 
normalized to the canonical value I = 174' /(47r). For 
comparison, we plot also (dashed line) the electric cur- 
rent distribution of a relativistic magnetic split monopole 
with the same amount of magnetic flux crossing the light 
cylinder (Michel 1991). 
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Fig. 5. Same as Fig. 4 with some amount of equatorial Fig. 7. The spindown torque parameter / (Eq. [T]) as a 
magnetospheric reconnection that corresponds to e = 0.4 function of the dissipation parameter e. In the absence of 
(Eq. [IDJ). reconnection (e = 0), / = 0.6 (CKF). 
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Fig. 6. Same as Fig. 4 with maximum equatorial magne- 
tospheric reconnection that corresponds to e = 1. 



